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ABSTRACT 

We give a review of different properties of dark matter halos. Taken from different 
publications, we present results on (1) the mass and velocity functions, (2) density 
and velocity profiles, and (3) concentration of halos. The results are not sensitive to 
parameters of cosmological models, but formally most of them were derived for popular 
flat ACDM model. In the range of radii r = (0.005 — l)r v ; r the density profile for 
a quiet isolated halo is very accurately approximated by a fit suggested by Moore et 
al.(1997): p oc 1/x {1 + x 15 ), where x = r/r s and r s is a characteristic radius. The 
fit suggested by Navarro et al.(1995) p oc l/x(l + x) 2 , also gives a very satisfactory 
approximation with relative errors of about 10% for radii not smaller than 1% of the 
virial radius. The mass function of z = halos with mass below ~ 1O 13 /i~ 1 M0 is 
approximated by a power law with slope a = —1.85. The slope increases with the 
redshift. The velocity function of halos with V maiX < 500 km/s is also a power law with 
the slope (3 = —3.8 — 4. The power-law extends to halos at least down to 10 km/s. It is 
also valid for halos inside larger virialized halos. The concentration of halos depends on 
mass (more massive halos are less concentrated) and environment, with isolated halos 
being less concentrated than halos of the same mass inside clusters. Halos have intrinsic 
scatter of concentration: at la level halos with the same mass have A(logc v i r ) = 0.18 
or, equivalently, AV max /V max = 0.12. Velocity anisotropyfor both subhalos and the 
dark matter is approximated by f3(r) = 0.15 + 2x/[x 2 + 4], where x is radius in units of 
the virial radius. 

Subject headings: cosmology:theory - galaxy structure - methods: numerical 



1. Introduction 

During the last decade there was an increasingly growing interest in testing predictions of 
variants of the cold dark matter (CDM) models at subgalactic (< 100 kpc) scales. This interest 
was first induced by indications that observed rotation curves in the central regions of dark matter 
dominated dwarf galaxies are at odds with predictions of hierarchical models. Specifically, it was 
argued (Moore 1994; Flores & Primack 1994) that circular velocities, v c (r) = [GM(< r)/r] 1 / 2 , 
at small galactocentric radii predicted by the models are too high and increase too rapidly with 
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increasing radius compared to the observed rotation curves. The steeper than expected rise of v c (r) 
implies that the shape of the predicted halo density distribution is incorrect and/or that the DM 
halos formed in CDM models are too concentrated (i.e., have too much of their mass concentrated 
in the inner regions). 

In addition to the density profiles, there is an alarming mismatch in predicted abundance of 
small-mass (< 10 8 — 10 9 /i _1 M©) galactic satellites and the observed number of satellites in the 
Local Group (Kauffmann, White & Guiderdoni 1993; Klypin et al. 1999; Moore et al. 1999). 
Although this discrepancy may well be due to feedback processes such as photoionization that 
prevent gas collapse and star formation in the majority of the small-mass satellites (e.g., Bullock, 
Kravtsov & Weinberg 2000), the mass scale at which the problem sets in is similar to the scale 
in the spectrum of primordial fluctuations that may be responsible for the problems with density 
profiles. In the age of precision cosmology that forthcoming MAP and Planck cosmic microwave 
background anisotropy satellite missions are expected to bring, tests of the cosmological models at 
small scales may prove to be the final frontier and the ultimate challenge to our understanding of 
the cosmology and structure formation in the Universe. However, this obviously requires detailed 
predictions and checks from the theoretical side and higher resolution/quality observations and 
thorough understanding of their implications and associated caveats from the observational side. 
In this paper we focus on the theoretical predictions of the density distribution of DM halos and 
some problems with comparing these predictions to observations. 

A systematic study of halo density profiles for a wide range of halo masses and cosmologies 
was done by Navarro, Frenk & White (1996, 1997; hereafter NFW), who argued that analytical 
profile of the form p(r) = p s (r /r s )~ l (\ + r/r s )~ 2 provides a good description of halo profiles in 
their simulations for all halo masses and in all cosmologies. Here, r s is the scale radius which, 
for this profile corresponds to the scale at which dlog p(r) /dlogr\ r=rg = —2. The parameters of 
the profile are determined by the halo's virial mass M v i r and concentration defined as c = r v i r /r s . 
NFW argued that there is a tight correlation between c and M v ; r , which implies that the density 
distributions of halos of different masses can in fact be described by a one-parameter family of 
analytical profiles. Further studies by Kravtsov, Klypin & Khokhlov (1997), Kravtsov et al. (1999), 
Jing (1999), Bullock et al. (2000), although confirming the c(M v ; r ) correlation, indicated that there 
is a significant scatter in the density profiles and concentrations for DM halos of a given mass. 

Following the initial studies by Moore (1994) and Flores k. Primack (1994), Kravtsov et al. 
(1999) presented a systematic comparison of the results of numerical simulations with rotation 
curves of a sample of seventeen dark matter dominated dwarf and low surface brightness (LSB) 
galaxies. Based on these comparisons, we argued that there does not seem to be a significant 
discrepancy in the shape of the density profiles at the scales probed by the numerical simulations 
(> 0.02 — 0.03r v i r , where r v ; r is halo's virial radius). However, these conclusions were subject to 
several caveats and had to be tested. First, observed galactic rotation curves had to be re-examined 
more carefully and with higher resolution. The fact that all of the observed rotation curves used 
in earlier analyses were obtained using relatively low-resolution HI observations, required checks 
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of the possible beam smearing effects. Also, a possibility of non-circular random motions in the 
central regions that could modify the rotation velocity of the gas (e.g., Binney & Tremain 1987, 
p. 198) had to be considered. Second, the theoretical predictions had to be tested for convergence 
and extended to scales < 0.01r v i r . 

Moore et al. (1998; see also a more recent convergence study by Ghigna et al. 1999) presented a 
convergence study and argued that mass resolution has a significant impact on the central density 
distribution of halos. They argued that at least several million particles per halo are required 
to model reliably the density profiles at scales < 0.01r v ; r . Based on these results, Moore et al. 
(1999) advocated a density profile of the form p(r) oc (r/ro)~ L5 [l + (r/Vo) 1 ' 5 ] -1 , that behaves 
similarly (p oc r~ 3 ) to the NFW profile at large radii, but is steeper at small r: p oc r~ L5 . Most 
recently, Jing & Suto (2000) presented a systematic study of density profiles for halo masses range 
2 x 10 12 /i _1 Mq — 5 x 1O 14 /i _1 M0. The study was uniform in mass and force resolution featuring 
~ 5 — 10 x 10 5 particles per halo and force resolution of « 0.004r v ; r . They found that galaxy-mass 
halos in their simulations are well fitted by profile 1 p(r) oc (r/ro)~ 1 ' 5 [l + r/ro]~ L5 , but that cluster- 
mass halos are well described by the NFW profile, with logarithmic slope of the density profiles at 
r = 0.01r vir changing from ss -1.5 for M vir ~ 10 12 /i _1 M o to « -1.1 for M vir ~5x lO 14 /t _1 M . 
Jing & Suto interpreted these results as evidence that profiles of DM halos are not universal. 

Rotation curves of a number of dwarf and LSB galaxies have recently been reconsidered using 
Ha observations (e.g., Swaters, Madore & Trewhella 2000; van den Bosch et al. 2000). The 
results show that for majority of galaxies Ha rotation curves are significantly different in their 
central regions than the rotation curves derived from HI observations. This indicates that the HI 
rotation curves are affected by beam smearing (Swaters et al. 2000). It is also possible that some 
of the difference may be due to real differences in kinematics of the two tracer gas components 
(ionized and neutral hydrogen). Preliminary comparisons of the new Ha rotation curves with 
model predictions show that NFW density profiles are consistent with the observed shapes of 
rotation curves (van den Bosch 2000; Navarro & Swaters 2000). Moreover, cuspy density profiles 
with inner logarithmic slopes as steep as ~ —1.5 also seem to be consistent with the data (van den 
Bosch 2000). Nevertheless, CDM halos appear to be too concentrated (Navarro k, Swaters 2000; 
McGaugh et al. 2000; Navarro & Steinmetz 2000) as compared to galactic halos, and therefore the 
problem remains. 

New observational and theoretical developments show that comparison of model predictions to 
the data is not straightforward. Decisive comparisons require reaching convergence of theoretical 
predictions and understanding the kinematics of the gas in the central regions observed galaxies. 
In this paper we present convergence tests designed to test effects of mass resolution on the density 
profiles of halos formed in the currently popular CDM model with cosmological constant (ACDM) 
and simulated using the multiple mass resolution version of the Adaptive Refinement Tree code 



1 Note that his profile is somewhat different than profile advocated by Moore et al., but behaves similarly to the 
latter at small radii. 
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(ART). We also discuss some caveats in drawing conclusions about the density profiles from the 
fits of analytical functions to numerical results and their comparisons to the data. 

2. Dark Matter Halos: the NFW and the Moore et al. profiles 

Before we fit the analytical profiles to real dark matter halos or compare them with observed 
rotational curves, it is instructive to compare different analytical approximations. Although the 
NFW and Moore et al. profiles predict different behavior of p(r) in the central regions of a halo, the 
scale where this difference becomes significant depends on the specific values of halo's characteristic 
density and radius. Table 2 presents different parameters and statistics associated with the two 
analytical profiles. For the NFW profile more information can be found in Klypin et al. (1999), 
Lokas & Mamon (2000), and in Widrow (2000). 

Each profile is set by two independent parameters. We choose these to be the characteristic 
density po an d radius r s . In this case all expressions describing different properties of the profiles 
have simple form and do not depend on concentration. The concentration or the virial mass appear 
only in the normalization of the expressions. The choice of the virial radius (e.g., Lokas & Mamon 
2000) as a scale unit results in more complicated expressions with explicit dependence on the 
concentration. In this case, one also has to be careful about definition of the virial radius, as there 
are several different definitions in the literature. For example, it is often defined as the radius, r2oo> 
within which the average density is 200 times the critical density. In this paper the virial radius is 
defined as radius within which the average density is equal to the density predicted by the top-hat 
model: it is 5th times the average matter density in the Universe. For the fio = 1 case the two 
existing definitions are equivalent. In the case of 0.$ = 0.3 models, however, the virial radius is 
about 30% larger than r2oo- 

There is no unique way of defining a consistent concentration for the different analytical profiles. 
Again, it is natural to use the characteristic radius r s to define the concentration: c = r v i r /r s . 
This simplifies expressions. At the same time, if we fit the same dark matter halo with the two 
profiles, we will get different concentrations because the values of corresponding r s will be different. 
Alternatively, if we choose to match the outer parts of the profiles (say, r > r s ) as closely as possible, 
we may choose to change the ratio of the characteristic radii r S) NF\v/ r s,Moorc in such a way that 
both profiles reach the maximum circular velocity v c i rc at the same physical radius r max . In this 
case, the formal concentration of the Moore et al. profile is 1.72 times smaller than that of the 
NFW profile. Indeed, with this normalization profiles look very similar in the outer parts as one 
finds in Figure 1. Table 2 also gives two other "concentrations". The concentration C1/5 is defined 
as the ratio of virial radius to the radius, which encompasses 1/5 of the virial mass (Avila- Reese et 
al. 1999). For halos with Cnfw ~ 5.5 this 1/5 mass concentration is equal to Cnfw- One can also 
define the concentration as the ratio of the virial radius to the radius at which the logarithmic slope 
of the density profile is equal -2. This scale corresponds to r s for the NFW profile and « 0.35r s for 
the Moore et al. profile. 
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Table 1. Comparison of NFW and Moore et al. profiles 



Parameter 


NFW Moore et al. 


Density 

x = r/r s 


P x(l + x) 2 

p oc x~ 3 for x S> 1 
p oc x _1 for x <C 1 
p/po = 1/4 at x = 1 
p/p = 1/21.3 at x = 2.15 


P ~ x Lb (l+x) Lh 
p oc x -3 for x ^ 1 

p oc x^ 1 ' 5 for x <C 1 

p/po = 1/2 at x = 1 

p/po = 1/3.35 at x = 1.25 


Mass 

M = 47rp r 3 J(x) 

= M vir /(x)//(C) 
M vir = f Pci n Q 6 th r* iT 


/(x) = ln(l + x)- T ^ 


f{x) = | ln(l + x 3 / 2 ) 


Concentration 
C = r v - ir / r s 


CnfW = l-72CMoore 

(for the same M^- 1T and r max ) 


CMoore = CnFw/1-72 
(~< , CMoore 


x / 5 ~ 0.86/ (C NFW ) + 0.1363 

(error <sx CMoore 
terror ^~ 3/10 

l°Moore ± J 

C 7= __2 = Cnfw 


[(l + < 0rC ) 1/5 -l] 2/3 

C-y=-2 = 2 3 / 2 CMoore 

~ 2.83CMoore 


Circular Velocity 

2 GM vil C f{x) 
Vcirc r vir xf(C) 

2 -^max fv^) 

x f(%ma,x) 

2 CM vir 
V ■ = 

vir . 

' vir 


^max ~ 2.15 

c 

<ax~0.216^ ir ^^ 


^max ~ 1.25 

c 

v 2 « 0.466-u 2 - 

Vax "^"""vir /(C) 
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Figure 1 presents comparison of the analytic profiles normalized to have the same virial mass 
and the same radius r max . We show results for halos of low and high values of concentration 
representative of cluster- and low-mass galaxy halos, respectively. The bottom panels show the 
profiles, while the top panels show corresponding logarithmic slope as a function of radius. The 
figure shows that the two profiles are very similar throughout the main body of the halos. Only in 
the very central region the differences become significant. The difference is more apparent in the 
logarithmic slope than in the actual density profiles. Moreover, for galaxy-mass halos the difference 
sets in at a rather small radius < 0.01r v i r , which would correspond to scales < 1 kpc for the typical 
dark matter dominated dwarf and LSB galaxies. In most analyses involving galaxy-size halos, the 
differences between NFW and Moore et al. profiles are irrelevant, and NFW profile should provide 
an accurate description of the density distribution. 

Note also that for galaxy-size (e.g., high-concentration) halos the logarithmic slope of the 
NFW profile does not reach its asymptotic inner value of —1 at scales as small as 0.01r v i r . For 
~ 10 12 h~ 1 M Q halos the logarithmic slope of the NFW profile is « —1.4 — 1.5, while for cluster-size 
halos this slope is ~ —1.2. This dependence of slope at a given fraction of the virial radius on 
the virial mass of the halo is very similar to the results plotted in Figure 3 of Jing & Suto (2000). 
These authors interpreted it as evidence that halo profiles are not universal. It is obvious, however, 
that their results are consistent with NFW profiles and the dependence of the slope on mass can 
be simply a manifestation of the well-studied c v i r (M) relation. 

To summarize, we find that the differences between the NFW and the Moore et al. profiles 
are very small (Ap/p < 10%) for radii above 1% of the virial radius. The differences are larger for 
halos with smaller concentrations. In case of the NFW profile, the asymptotic value of the central 
slope 7 = — 1 is not achieved even at radii as small as l%-2% of the virial radius. 

3. Properties of Dark Matter Halos 

Some properties of halos depend on large-scale environment in which the halos are found. We 
will call a halo distinct if it is not inside a virial radius of another (larger) halo. Halo is called 
subhalo, if it is inside of another halo. The number of subhalos depends on the mass resolution 
- the deeper we go, the more subhalos we will find. Most of the results below are based on a 
simulation, which was complete to masses down to 1O 11 /i _1 M or, equivalently, to the maximum 
circular velocity of 100 km/s. 

Mass and Velocity Distribution functions. Extensive analysis of halo mass and velocity 
function was done by Sigad et al. (2000) for halos in the ACDM model. Additional results can 
also be found in Ghigna et al. (1999); Moore et al. (1999b); Klypin et al. (1999); Gottlober, 
Klypin & Kravtsov (1998). Figure 2 compares the mass function of subhalos and distinct halos. 
The Press-Schechter approximation overestimates the the mass function by a factor of 2 for M < 
5 x 1O 12 /i _1 M and it somewhat underestimates it at larger masses. More advanced approximation 
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Fig. 1. — Comparison of the Moore et al. and the NFW profiles. Each profile is normalized to have the 
same virial mass and the same radius of the maximum circular velocity. Left panels: High concentration 
halo with concentration typical for small galaxies Cnfw — 17 Right panels: Low concentration halo with 
concentration typical for clusters of galaxies. The deviations are very small (< 3%) for radii r > l/2r s . Top 
panels show the local logarithmic slope of the profiles. Note that for the high concentration halo the slope 
of the profile is significantly larger than the asymptotic value -1 even at very small radii r rs 0.01/r V j r . 

given by Sheth & Tormen is more accurate. On scales below 1O 14 /i _1 M0 the mass function is close 
to a power law with the slope a ~ —1.8. There is no visible difference in the slope for subhalos and 
for the distinct halos. 

For each halo one can measure the maximum circular velocity Vmax- In many cases (especially 
for subhalos) V max is a better measure of how large is the halo. It is also closer related with 
observed properties of galaxies hosted by halos. Figure 3 presents the velocity distribution function 
of different types of halos. In addition to distinct halos and subhalos, we show also isolated halos 
and halos in groups and clusters. Here isolated halos are defined as halos with mass less than 
10 13 fc -1 M©, which are not inside a larger halo and which do not have subhalos more massive than 
10 u h^ 1 M & . The velocity function is approximated by a power law dn = $*V^axdVWx with the 
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Fig. 2. — Mass function for distinct halos (top) and for subhalos bottom). Raw counts are marked by 
symbols with error bars. The curves are Schechter-function fits. Press-Schechter (dotted) and Sheth-Tormen 
(dashed) predictions for distinct halos are also shown. On scales below 10 14 h~ 1 M ( ? ) the mass function is 
close to a power law with the slope a « —1.8. There is no visible difference in the slope for subhalos and for 
the distinct halos. (After Sigad et al. (2000)) 



slope (3 ~ —3.8 for distinct halos. The slope depends on environment: (3 ~ —3.1 for halos in groups 
and (3 —4 for isolated halos. Klypin et al. (1999) and Ghigna et al. (1999) found that the 
slope j3 ~ —3.8 — 4 of the velocity function extends to much smaller halos with velocities down to 
10 km/s. 

Correlation between characteristic density and radius. The halo density profiles arc 
approximated by the Navarro-Frenk- White profile: 



P (r/r )[l + r/r } 2 U 

Kravtsov et al. (1999) find the correlation between the two parameters of halos po and r s . 
Fifure 4 compares results for the dark matter halos with those for the the dark matter dominated 
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Fig. 3. — Velocity functions for isolated halos (squares) and for halos in groups and clusters. Halos with 
mass less than 10 13 /i _1 M© are used for the plots. (After Sigad et al. (2000)) 

Low Surface Brightness (LSB) galaxies and dwarf galaxies. Halos are consistent with observational 
data: smaller halos are denser. 

Correlations between mass, concentration, and redshift. Navarro et al. (1997) argued 
that the halo profiles have a universal shape in the sense that profile is uniquely defined by virial 
mass of the halo. Bullock et al. (2000) analyzed concentrations of thousands of halos at different 
redshifts. To some degree they confirm conclusions of Navarro et al. (1997): halo concentration 
correlates with its mass. But some significant deviations were also found. There is no one-to- 
one relation between concentration and mass. It appears the the universal profile should only 
be treated as a trend: halo concentration does increase as the halo mass decreases, but there 
are large deviations for individual halos from that "universal" shape. Halos have intrinsic scatter 
of concentration: at la level halos with the same mass have A(logc v ; r ) = 0.18 or, equivalently, 

A^max/^max — 0.12. 

Velocity anisotropy. Inside a large halo subhalos or dark matter particles do not move on 
either circular or radial orbits. Velocity ellipsoid can be measured at each position inside halo. 
It can be characterized by anisotropy parameter defined as /3(r) = 1 — V^/2V^. Here V"j_ is the 
velocity dispersion perpendicular to the radial direction and V? is the radial velocity dispersion. 
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Fig. 4. — Correlation of the characteristic density po and radius ro for the dwarf and LSB galaxies (solid and 
open circles) and for the dark matter halos (crosses) in different cosmological models. Halos are consistent 
with observational data: smaller halos are denser. (After Kravtsov et al. (1999)) 

For pure radial motions (3 = 1. For isotropic velocities (3 = 0. Function (3{r) was estimated for 
halos in different cosmological models ( see Colin et al. (1999) for references). By studying 12 rich 
clusters with many subhalos inside each of them, Colin et al. (1999) found that both the subhalos 
and dark matter particles can be described by the same anisotropy parameter 

2x 

(3{r) =0.15 + ^^, x = r/r viT (2) 
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Fig. 5. — Left Panel Dependence of concentration with mass for distinct halos. The thick solid curve is the 
median value. The errors are errors of the mean due to sampling. The outer dot-dashed curves encompass 
68% of halos in the simulations. The dashed curves and arrows indicate values corrected for the noise in halo 
profiles. Thin curves are different analytical models. Right Panel Median halo concentration as a function of 
mass for different redshifts. The thin lines show predictions of an analytical model. (After Bullock et al. 
(2000)). 

4. Halo profiles: Convergence study 

The following results are based on Klypin et al. (2000) 

4.1. Numerical simulations 

Using the ART code (Kravtsov, Klypin & Khokhlov 1997; Kravtsov 1999), we simulate a flat 
low-density cosmological model (ACDM) with = 1 — = 0.3, the Hubble parameter (in units 
of 100 kms _1 Mpc _1 ) h = 0.7, and the spectrum normalization as = 0.9. We have run two sets of 
simulations with 30/i _1 Mpc and 25/i _1 Mpc computational box. The first simulations were run to 
the present moment z = 0. The second set of simulations had higher mass resolution and therefore 
produced more halos but were run only to z = 1. 

In all of our simulations step in the expansion parameter was chosen to be Aan = 2 x 10 -3 on the 
zero level of resolution. This gives about 500 steps for an entire run to z = 0. A test run was done 
with twice smaller time-step for a halos of comparable mass (but with smaller number of particles) 
as studied in this paper. We did not find any visible deviations in the halo profile. In the first set of 
simulations, the highest level of refinement was ten, which corresponds to 500 x 2 10 ~ 500, 000 time 
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steps at the tenth level. For the second set of simulation, nine levels of refinement were reached 
which corresponds to 128, 000 steps at the ninth level. 

In the following sections we present results for four halos. The first halo (A) was the only halo 
selected for resimulation in the first set of simulations. In this case the selected halo was relatively 
quiescent at z = and had no massive neighbors. The halo was located in a long filament bordering 
a large void. It was about 10 Mpc away from nearest cluster-size halo. After the high-resolution 
simulation was completed we found that the nearest galaxy-size halo was about 5 Mpc away. The 
halo had a fairly typical merging history with M(t) track slightly lower than the average mass 
growth predicted using extended Press-Schechter model. The last major merger event occured at 
z ~ 2.5; at lower redshifts the mass growth (the mass in this time interval has grown by a factor of 
three) was due to slow and steady mass accretion. 

The second set of simulations was done in a different way. In the low resolution run we selected 
three halos in a well pronounced filament. Two of the halos are neighbors located at about 0.5 Mpc 
from each other. The third halo was 2 Mpc away from this pair. Thus, the halos were not selected 
to be too isolated as was the case in the first set of runs. Moreover, the simulation was analyzed at 
an earlier moment (z = 1) where halos are more likely to be unrelaxed. Therefore, we consider the 
halo A from the first set as an example of a rather isolated well relaxed halo. In many respects, 
this halo is similar to halos simulated by other research groups that used multiple mass resolution 
techniques. The three halos from the second set of simulations can be viewed as representative of 
more typical halos, not necessarily well relaxed and located in more crowded environments. 

Parameters of the simulated dark matter halos are listed in Table 2. Columns in the table 
present (1) Halo "name" (halos Ai, A2, A3 are the halo A re-simulated times with different res- 
olutions); (2) redshift at which halo was analyzed; (3-5) virial mass, comoving virial radius, and 
maximum circular velocity. At z = (z = 1) the virial radius was estimated as the radius within 
which the average overdensity of matter is 340 (180) times larger than the mean cosmological den- 
sity of matter at that redshift; (6) the number of particles within the virial radius; (7) the smallest 
particle mass in the simulation; (8) formal force resolution achieved in the simulation. As we will 
show below, convergent results are expected at scales larger than four times the formal resolution; 
(9) halo concentration as estimated from NFW profile fits to halo density profiles; (10) maximum 
relative error of the NFW fit: /ONFw/Phalo — 1 (the error was estimated inside 50/i _1 kpc radius); 
(11) the same as in the previous column, but for the fits of profile advocated by Moore et al. 

Halo A in the first set of simulations was re-simulated three times with increasing mass resolu- 
tion. For each simulation, we considered outputs at four time moments in the interval to z = 0—0.03. 
Parameters of the halos in these simulations averaged over the four moments are presented in the 
first three rows of the Table 2. We do not find any systematic change with resolution in the values 
of halo parameters both on the virial radius scale and around the maximum of the circular velocity 
(r = (30-A0)h- l kpc). 

Left panel in Figure 6 shows the central region of the halo Ai (see Table 2). This plot is similar 
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to the Fig. la in Moore et al. (1998) in that all profiles are drawn to the formal force resolution. The 
straight lines indicate slopes of two power-laws: 7 = — 1 and 7 = —1.4. The figure indeed shows 
that at around 1% of the virial radius the slope is steeper than -1 and the central slope increases 
as we increase the mass resolution. Moore et al. (1998) interpreted this behavior as evidence that 
profiles are steeper than predicted by the NFW profile. We also note that the results of our highest 
resolution run A\ are qualitatively consistent with results of Kravtsov et al. (1998). Indeed, if the 
profiles are considered down to the scale of two formal resolutions, the density profile slope in the 
very central part of the profile r < 0.01r v i r is close to 7 = —0.5. 

The profiles in Figure 6 reflect the density distribution in the cores of simulated halos. However, 
the interpretation of these profiles is not straightforward because it requires assessment of numerical 
effects. The formal resolution usually does not even correspond to the scale where the numerical 
force is fully Newtonian (usually it is still considerably "softer" than the Newtonian value). In the 
ART code, the interparticle force reaches (on average) the Newtonian value at about two formal 
resolutions (see Kravtsov et al. 1997). The effects of force resolution can be studied by resimulating 
the same objects with higher force resolution and comparing the density profiles. Such convergence 
study was done in Kravtsov et al. (1998) where it was found that for a fixed mass resolution halo 
density profiles converge at scales above two formal resolutions. Second, local dynamical time for 
particles moving in the core of a halo is very short. For example, particles on the circular orbit of 
the radius l/i _1 kpc from the center of halo A makes about 200 revolutions over the Hubble time. 
Therefore, if the time step is insufficiently small, numerical errors in these regions will tend to 
grow especially fast. The third possible source of numerical errors is mass resolution. Poor mass 
resolution in simulations with good force resolution may, for example, lead to two-body effects (e.g., 
Knebe et al. 2000). Insufficient number of particles may also result in "grainy" potential in halo 
cores and thereby affect accuracy of orbit integration. In these effects, the mass resolution may be 
closely inter-related with force resolution. 

It is clear thus that in order to make conclusions not affected by numerical errors, one has 
to determine the range of trustworthy scales using convergence analysis. Right panel in Figure 6 
shows that for the halo A simulations the convergence for vastly different mass and force resolution 
is reached for scales > 4 formal force resolutions (all profiles in this figure are plotted down to the 
radius of 4 formal force resolutions). For all resolutions, there are more than 200 particles within 
the radius of four resolutions from the halo center. For the highest resolution simulation (halo Ai) 
the convergence is reached at scales > 0.005r v ; r . 

In order to judge which profile provides a better description of the simulated profiles we fitted 
the NFW and Moore et al. analytic profiles. Figure 7 presents results of the fits and shows that 
both profiles fit the numerical profile equally well: fractional deviations of the fitted profiles from 
the numerical one are smaller than 20% over almost three decades in radius. It is clear thus that 
the fact that numerical profile has slope steeper than —1 at the scale of ~ 0.01r v j r does not mean 
that good fit of the NFW profile (or even analytic profiles with shallower asymptotic slopes) cannot 
be obtained. 
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Fig. 6. — Left panel Density profiles of halo A simulated with different mass and force resolutions. The 
profiles are plotted down to the formal force resolution of each simulation. Right panel The profiles plotted 
down to four formal resolutions. It is clear that for vastly different mass (from 2000 to 120000 particles in 
the halo) and force (from 3.66/i _1 kpc to 0.23/i _1 kpc) resolutions the convergence is reached at these scales. 

There is certainly a certain degree of degeneracy in fitting various analytic profile to numerical 
results. Figure 8 illustrates this further by showing results of fitting profiles (solid lines) of the form 
p(r) oc (r/r )- 7 [l + (r/r ) a ]- (/3 - a)/7 to the same (halo A±) simulated halo profile shown as solid 
circles. The legend in each panel indicates the corresponding values of a, (3, and 7 of the fit; the 
digit in parenthesis indicates whether the parameter was kept fixed (0) or not (1) during the fit. 
The right two panels show fits of the NFW and Moore et al. profiles; the bottom left panel shows 
fit of the profiles used by Jing & Suto (2000). The top left panel shows a fit in which the inner 
slope was fixed but a and (3 were fit. The figure shows that all four analytic profiles can provide a 
nice fit to the numerical profile in the whole range 0.005 — lr v ; r . 



4.2. Halo profiles at z = 1 

As we mentioned in § 4.1, the halo A analyzed in the previous section is somewhat special 
because it was selected as an isolated relaxed halo. In order to reach unbiased conclusions, in this 
section we will present analysis of halos from the second set of simulations (halos B, C, and D in 
Table 2) which were not selected to be relaxed or isolated. Based on the results of convergence study 
presented in the previous section, we will consider profiles of these halos only at scales above four 
formal resolutions use results starting only from 4 formal resolutions and not less than 200 particles. 
Note that these conditions are probably more stringent than necessary because these halos were 
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Fig. 7. — Fits of the NFW and Moore et al. halo profiles to the profile of halo Ai ) bottom panel. The 
top panel shows fractional deviations of the analytic fits from the numerical profile. Note that both analytic 
profiles fit numerical profile equally well: fractional deviations are smaller than 20% over almost three decades 
in radius. 

simulated with 5 — 7 times more particles per halo. There is an advantage in analyzing halos at 
a relatively high redshift. Halos of a given mass will have lower concentration (see Bullock et al. 
2000). Lower concentration implies a large scale at which the asymptotic inner slope is reached. 
Profiles of the high redshift halos should therefore be more useful in discriminating between the 
analytic models with different inner slopes. 

We found that substantial substructure is present inside the virial radius in all three halos. 
Figure 9 shows profiles of these halos at z = 1. There profiles are not as smooth as that of halo 
Ai due to the substructure. Note that bumps and depressions visible in the profiles cannot are 
significantly larger amplitude than the shot noise. Halo C appeared to be the most relaxed of the 
three halos. This halo had the last major merger somewhat earlier than the other two. Halo D 
had a major merger event at z ~ 2. Remnant of the merger are still visible as a hump at radii 
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Fig. 8. — Analytic fits to the density profile of the halo A± (see Table 2) from our set of simulations. The fits 
are of the form p(r) oc (r/ro) '[1 + (r/ro)"] - '-' 3- "^ 7 . The legend in each panel indicates the corresponding 
values of a, (3, and 7 of the fit; the digit in parenthesis indicates whether the parameter was kept fixed (0) or 
not (1) during the fit. Note that various sets of parameters a, (3, 7 provide equally good fits to the simulated 
halo profile in the whole range resolved range of scales w 0.005 — lr v ; r . This indicates a large degree of 
degeneracy in parameters a, (3, and 7 

around 100/i _1 kpc. The non-uniformities of profiles cause by substructure may substantially bias 
analytic fits to the entire range of scales below the virial radius. Therefore, we used only the central, 
presumably more relaxed, regions in the analytic fits: r < 50/i _1 kpc for halo D and r < 100/i -1 kpc 
for halos B and C (fits using only central 50/i _1 kpc did not change results). 
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The best fit parameters were obtained by minimizing the maximum fractional deviation of 
the fit: max(abs(log/9g t ) — logphalo)- Minimizing the sum of squares of deviations (x 2 )> as is often 
done, can result in larger errors at small radii with the false impression that the fit fails because it 
has a wrong central slope. The fit that minimizes maximum deviations, improves the NFW fit for 
points in the range of radii (5 — 20)/i _1 kpc, where the NFW fit would appear to be below the data 
points if the fit was done by the \ 2 minimization. This improvement comes at the expense of few 
points around l/i -1 kpc. For example, if we fit halo B by using x 2 minimization, the concentration 
decreases from 12.3 (see Table 2) to 11.8 We also made a fit for halo B assuming even more stringent 
limits on the effects of numerical resolution. By minimizing the maximum deviation we fitted the 
halo starting at six times the formal resolution. Inside this radius there were about 900 particles. 
Resulting parameters of the fit were close to those in Table 2: Cnfw = H-8, and maximum error 
of the NFW fit was 17%. 

We found that the errors in the Moore et al. fits were systematically smaller than those of the 
NFW fits, though the differences were not dramatic. Moore et al. fit failed for halo D. It formally 
gave very small errors, but this was done for a fit with unreasonably small concentration C = 2. 
When we constrained the approximation to have about twice larger concentration as compared 
with the best NFW fit, we were able to obtain a reasonable fit (this fit is shown in Figure 9). 
Nevertheless, the central part was fit poorly in this case. 

Our analysis therefore failed to determine which analytic profile provides a better description 
of the density distribution in simulated halos. Despite the larger number of particles per halo and 
lower concentrations of halos, results are still inconclusive. Moore at al. profile is a better fit to the 
profile of halo C; the NFW profile is a better fit to the central part of halo D. Halo B represents an 
intermediate case where both profiles provide equally good fits (similar to the analysis of halo A) . 

Note that there seem to be real deviations in parameters of halos of the same mass. Halos B 
and D have the same virial radii and nearly the same circular velocities, yet their concentrations 
are different by 30%. We find the same differences in estimates of C1/5 concentrations, which do 
not depend on specifics of an analytic fit. The central slope at around lkpc also changes from halo 
to halo. 

5. Summary 

We run a series of simulations with vastly different mass and force resolution with the goal of 
studying the shape of the dark matter halo profile in central parts of galaxy-size halos. We use 
a modified version of the ART code, which is capable of handling particles with different masses, 
variable force, and time resolution. In runs with the highest resolution, the code achieved (formal) 
dynamical range of 2 17 = 131,072 with 500,000 steps for particles at highest level of resolution. 

Our conclusions regarding the convergence of the profiles differ from those of Moore et al. 
(1997). If we take into account only radii, at which we believe numerical effects (the force resolution, 
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Fig. 9. — Profiles of halos B, C, D at z = 1. Profiles of halos C and D were offset downwards by factors 
of 10 and 100 for clarity. Solid curves show simulated profiles, while dotted and dot-dashed curves show the 
NFW and Moore et al. fits, respectively. The halo profiles in simulations are plotted down to four formal 
resolution. Each halo had more than 200 particles inside the smallest plotted scale. 

the resolution of initial perturbations, and the two-body scattering) are small, then we find that the 
slope and the amplitude of the density do not change when we change the force and mass resolution. 
This result is consistent with what was found in simulations of the "Santa Barbara" cluster (Frenk 
et al., 1999): at a fixed resolved scale results do not change as the resolution increases. For the 
ART code the results converged at 4 times the formal force resolution and more than 200 particles. 
These limits of convergence very likely depend on particular code used and on the duration of 
integration. 

We reproduce results of Moore et al. regarding the convergence and results of Kravtsov et 
al. (1998) regarding shallow central profiles, but only when we consider points inside unresolved 
scales. We conclude that those results followed from overly optimistic interpretation of numerical 
accuracy of simulations. 
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For the galaxy-size halos considered in this paper with masses M v - ir = 7 x 1O 11 /i" 1 M — 2 x 
1O 12 /i _1 M0 and concentrations C = 9 — 17 both the NFW profile p oc r _1 (l +r)~ 2 and the Moore 
et al. profile p oc r~ L5 (l + r 1 - 5 ) -1 give good fits with accuracy about 10% for radii not smaller 
than 1% of the virial radius. None of the profiles is significantly better then the other. 

Halos with the same mass may have different profiles. No matter what profile is used - NFW 
or Moore et al. - there is no universal profile: just halo mass does not yet define the density profile. 
Nevertheless, the universal profile is extremely useful notion which should be interpreted as the 
general trend C(M) of halos with larger mass to have lower concentration. Deviations from the 
general C(M) are real and significant (Bullock et al. 2000). It is not yet clear, but seems very likely 
that the central slopes of halos also have real fluctuations. The fluctuations in the concentration 
and the central slopes are important for interpretation of the central parts of rotation curves. 

We acknowledge the support of the grants NAG- 5- 3842 and NST- 9802787. A.V.K. ac- 
knowledges support by NASA through Hubble Fellowship grant HF-01121.01-99A from the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. Computer simulations presented in this paper 
were done at the National Center for Supercomputing Applications (NCSA), Urbana-Champaign, 
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Table 2. Parameters of Halos 





z 


M vir 




^riiax 




m pa,rt 


Form. res. 


Cnfw 


RelEr 


RelEr 






M Q /h 


kpc/h 


km/s 




M Q /h 


kpc/h 




NFW 


Moore 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Ai 





1.97- 10 12 


257 


247.0 


1.2 • 10 5 


1.6 • 10 7 


0.23 


17.4 


0.17 


0.20 


A 2 





2.05 • 10 12 


261 


248.5 


1.5 • 10 4 


1.3 • 10 8 


0.91 


16.0 


0.13 


0.16 


A 3 





1.98 • 10 12 


256 


250.5 


1.9 • 10 3 


1.1 • 10 9 


3.66 


16.6 


0.16 


0.10 


B 


1 


8.5 • 10 11 


241 


195.4 


7.1 • 10 5 


1.2 • 10 6 


0.19 


12.3 


0.23 


0.16 


C 


1 


6.8 • 10 11 


208 


165.7 


5.0 • 10 5 


1.2 • 10 6 


0.19 


11.9 


0.37 


0.20 


D 


1 


9.6 • 10 11 


245 


202.4 


7.9 • 10 5 


1.2 • 10 6 


0.19 


9.5 


0.25 


0.60 



